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Abstract. Seventeen accreting neutron star pulsars, mostly high mass X-ray binaries with half of 
them Be-type transients, are known to exhibit Cyclotron Resonance Scattering Features (CRSFs) 
in their X-ray spectra, with characteristic line energies from 10 to 60keV. To date about two thirds 
of them, plus a few similar systems without known CRSFs, have been observed with Suzaku. We 
present an overview of results from these observations, including the discovery of a CRSF in the 
transient 1A 1118—61 and pulse phase resolved spectroscopy of GX 301—2. These observations 
allow for the determination of cyclotron line parameters to an unprecedented degree of accuracy 
within a moderate amount of observing time. This is important since these parameters vary - e.g., 
with orbital phase, pulse phase, or luminosity - depending on the geometry of the magnetic field 
of the pulsar and the properties of the accretion column at the magnetic poles. We briefly introduce 
a spectral model for CRSFs that is currently being developed and that for the first time is based 
on these physical properties. In addition to cyclotron line measurements, selected highlights from 
the Suzaku analyses include dip and flare studies, e.g., of 4U 1907+09 and Vela X-l, which show 
clumpy wind effects (like partial absorption and/or a decrease in the mass accretion rate supplied by 
the wind) and may also display magnetospheric gating effects. 
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netic fields — stars: neutron — X-rays: binaries 
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FIGURE 1. Left: Swift -BAT 15-50 keV light curve of the transient accreting pulsar 1A 11 18-61 
showing the 2009 outburst of the source. Vertical lines indicate the times of the January 15, 2009, and 
January 28, 2009, Suzaku pointings (exposures of ^50 ks and ^30 ks, respectively). Right: RXTE-ASM 
1-12 keV orbital profile of the persistent accreting pulsar GX301— 2. Vertical lines indicate the orbital 
phases of the August 25, 2008, and January 5, 2009, Suzaku pointings (exposures of ^10ks and ^60 ks, 
respectively). After [1] and [2]. 



1. INTRODUCTION 

Accreting X-ray pulsars were first discovered in the 1970s with the detection of the 4.8 s 
pulsation of Cen X-3 [3]. In these sources the compact object is a highly magnetized 
neutron star leading to accretion of material along the magnetic field lines and the for- 
mation of X-ray emitting accretion columns above the neutron star's magnetic poles. 
These sources can also show absorption-line-like scattering features in their X-ray spec- 
tra at multiples of E cyc i ~ 11.6keV x 5[10 12 G] since for such strong magnetic fields 
electron energies are quantized. Observations of sources with cyclotron resonant scat- 
tering features (CRSF, cyclotron lines) provide us with the most direct information for 
studying strong magnetic fields. Suzaku is especially well suited to find/study cyclotron 
lines due to its broadband X-ray coverage at high sensitivity with comparatively low 
background. With O- or B-type donor stars most CRSF sources are high mass X-ray 
binaries (exceptions: Her X-l, 4U 1626—67). Those in the first block of Table 1 (plus 
persistent X Per and some of the potential CRSF candidates in the third block) are Be 
type X-ray binaries which are generally transient [4]. They typically have eccentric or- 
bits and their X-ray activity is triggered when the neutron star crosses a disk of material 
that surrounds the equator of the Be/Oe star. Between outbursts they can show years long 
quiescence periods. While the CRSF sources listed in the second block of Table 1 are 
persistent X-ray emitters they are very variable on all observed time scales, reflecting 
stellar wind, magnetosphere (see, e.g., [5, 6] for Vela X-l), or accretion disk structure 
(see, e.g., [7] for Her X-l). Fig. 1 shows examples of typical long term variability for a 
transient and a persistent source. For recent reviews on accreting pulsars see [8] and [9]. 



TABLE 1. Properties of a sample of accreting pulsars. 1-9: All known transient CRSF sources. 10-17: 
All known persistent CRSF sources. 18-27: Selected potential CRSF candidates. All sources with names 
not printed in italics have been observed with Suzaku (1 x — 4x) at the time of writing. The pulse period 
'pulse usually varies with luminosity and time. The same is often true for the cyclotron line energies E cyc \ 
which can also vary with pulse phase. Here we quote representative values from the overall literature. 
Suzaku results for sources with names in bold face are discussed in this work. 
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FIGURE 2. Left: The main panel shows the Suzaku broad band spectrum and best fit of the brighter 
1A 1118—61 observation with an inset highlighting the iron line region. Residuals are shown including 
no CRSF (a), including one CRSF at 58keV (b), and including a second CRSF at 112keV (c, best fit). 
Right:. The main panel shows the Suzaku broad band spectrum and best fit of the brighter GX301— 2 
observation with an inset highlighting the emission line residuals from neutral material which remain 
after modeling Fe Ka. (The back illuminated XIS1 was modeled separately here, with consistent results.) 
Residuals are shown including no CRSF and not accounting for partial covering absorption (a), including 
a CRSF at 35 keV (b), and also accounting for partial covering (c, best fit). After [1] and [2]. 



2. SUZAKU OBSERVATIONS 
2.1. Discovery of a CRSF in 1A 1118+61 

After only two previously observed outbursts (1974, 1992) the Be X-ray binary 
1 A 1 1 18 — 61 went into outburst in January 2009 and was observed with Suzaku at 
the ~ 500mCrab outburst peak and two weeks later at the end of the main outburst 
(Fig. 1, left). Together with quasi-simultaneous RXTE observations [6] these observa- 
tions enabled the discovery of a ~ 55keV cyclotron line in this source [1] (Fig. 2, left). 
RXTE, Swift-BAT, and Suzaku also co-discovered another cyclotron line, namely the 
~ 54keV feature in GX304— 1 [25, 13]. Furthermore there are indications that the line 
in 1A 1118 — 61 may be luminosity dependent. With E cyc \ values of 58.2+^keV and 
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FIGURE 3. Left: Results from a - 200 ks RXTE observation of GX 301-2, after [27]. The upper 
panel shows the double peaked pulse profile in the 19-60 keV energy band (twice for clarity). The lower 
panel shows that the cyclotron line energy varies by ^20% over one cycle. Vertical lines indicate pulse 
profile minima. Right: Results from a ~ 60 ks Suzaku observation of GX 301—2 ("obs 2" from Fig. 1 and 
Fig. 2), after [2]. The red data points confirm the E cyc \ variations, translated into directly E cyc \ proportional 
B-field variations. The black curve shows the best fit magnetic dipole model described in the text. The blue 
histogram indicates the l-10keV pulse profile (arbitrary normalization). 

47.4+2 3 keV for the flare and decline observation, respectively, the source would show 
a positive L-E cyc \ correlation like Her X-l. Negative or no correlations have recently 
been observed for other CRSF sources [8]. A possible explanation for the existence of 
positive and negative correlations is the different state of the accretion column in the 
sub- and super- Eddington regimes [26]. 



2.2. Pulse Phase Resolved Analysis of GX 301—2 

The persistent accreting pulsar GX301— 2 was observed twice with Suzaku, during 
less bright and less well studied orbital phases outside of the pre-periastron flare (Fig. 1, 
right). These observations allowed for a detailed study of the highly absorbed phase 
averaged spectrum, including several fluorescence emission lines and the well known 
37 keV CRSF (Fig. 2, right) as well as of pulse phase resolved spectra for the second, 
longer observation [2]. Strong variations of the cyclotron line energy iicyci with pulse 
phase had been discovered in an RXTE observation dominated by the comparatively 
bright pre-periastron flare [27] (Fig. 3, left, lower panel). Several other CRSF sources 
are known to show similarly strong variations (see, e.g., [28] and references given by 
[27, 1]). Even though the Suzaku observation does not cover the bright pre-periastron 
flare a clear pattern for the ^cyci variations emerges within a comparatively modest 
exposure time (Fig. 3, right). Assuming a simple magnetic dipole the 5-field projected 
onto the line of sight can be calculated for each phase bin, depending on a set of 
geometric angles. This leads to four symmetrically equivalent best fit solutions. Model 
constraints derived using the observed ratio of the width of the cyclotron line and £cyci 
favor a solution wherein the spin axis is tilted by 15° with respect to the line of sight [2]. 
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FIGURE 4. Left: Energy-resolved, background-subtracted light curves (upper panel), and softness 
ratios (lower panel) from XISO for the 2006 Suzaku observation of 4U 1907+09. Right: Color-color 
diagram from added XISO and XIS3 data for the same observation. Boxed data points correspond to the 
deep dip in the light curve between 60 ks and 100 ks. Grid lines indicate hardness values from spectral 
models obtained by increasing the covering fraction at constant A^h (blue) or vice versa (red) compared to 
the best fit spectral model for times without additional partial covering absorption or flaring. From [14]. 

2.3. Dips & Flares: Stellar Wind Structure 

2.3.1. 4U 1907+09 

The ~20mCrab accreting pulsar 4U 1907+09 is characterized by a highly variable 
persistent flux as well, but in this case the dips and flares do not seem to be related to 
the orbital phase [29, 14]. Suzaku observed the source on May 2-3, 2006, and April 19- 
20, 2007, for respectively ~60ks and ~ 80 ks. Both observations display dips and flares 
(Fig. 4, left, shows the 2006 case). It had been noted before that the dipping might not be 
due to absorption but due to cessation of accretion [29]. Using color-color diagrams [14] 
found that while some of the flux and hardness variability is consistent with absorption 
there are indeed events where this is not the case, i.e., for the deep dip in the 2006 dataset 
(Fig. 4, right). They conclude that the overall behavior is consistent with a clumpy wind. 



2.3.2. VelaX-1 

Even more extreme flux variability, so called off states and giant flares - both with 
flux changes of up to a factor of ~20 - have been discovered with INTEGRAL for Vela 
X-l [30]. During the off states no pulsations were observed with INTEGRAL. Suzaku 
detected three off-sates during a 100 ks observation on June 17-18, 2008, lasting up 
to 3.2 ks [6]. Due to Suzaku's sensitivity pulsations could still be detected. [6] suggest 
that while the off states might be triggered by stellar wind inhomogeneities they are 
consistent with gated accretion where Kelvin-Helmholz instabilities allow for some 
matter to leak into the magneto sphere. 
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FIGURE 5. Comparison of spectra resulting from the previous (black) and current (red) version of the 
MC code for simulating physical CRSFs, for three different photon incident angles. From [31]. The input 
X-ray continuum spectrum is an empirical one (a power law with a Fermi-Dirac cutoff, characterized by 
EfoU — 10 keV and E cut = 12 keV). The lines, however, are produced in a cylindrical accretion column with 
a magnetic field of B/B C rit = 0.05 and an electron plasma temperature and optical depth of kT e — 3 keV 
and T es = 10~ 3 , respectively. Note that B crit = mlc 3 /efr = 4.414 x 10 13 G. 



3. A PHYSICALLY MOTIVATED CRSF MODEL 



Currently empirical models are used to describe cyclotron line features, parametrized, 
e.g., as absorption with a Gaussian optical depth profile or by a Lorentzian profile (gabs 
and cyclabs in xspec, respectively). Work has been going on since the 1990s to 
describe CRSFs in terms of physical parameters like the temperature kT e and optical 
depth T es of the electron plasma in the line forming region, as well as the 5-field strength 
and the accretion geometry [32, 33]. This endeavor reached a recent point of culmination 
with the implementation of a Monte Carlo (MC) code describing the scattering process, 
the calculation of a grid of Green's functions which can be used to apply the process to 
any continuum spectrum, and the development of the xspec model cyclomc for this 
task by [34]. First cyclomc applications to observed X-ray spectra, e.g., of V 0332+53, 
resulted in generally acceptable fits and parameter values, with the caveat of stronger 
emission wings than observed and different relative line depths than expected [34]. With 
a correction of the scattering cross sections in the most recent iteration of the MC code 
by [31] this has been resolved. Fig. 5 shows a comparison of CRSFs calculated with the 
two MC versions for three different photon incident angles with respect to the magnetic 
field in a cylindrical accretion column. The team is currently in the process of calculating 
a new grid of Green's functions to be used with cyclomc. 



4. SUMMARY AND CONCLUSIONS 



So far Suzaku has observed eleven of the seventeen cyclotron line sources at least 
once. Two were identified as CRSF sources for the first time (1A 1118 — 61, GX 304—1). 
Detailed phase resolved studies can be performed that allow us to constrain the 5-field 
geometry (e.g., GX 301-2). The important E cyc \ — L relationship can be studied (e.g., 
1A 1118—61, 1 A 0535+26), potentially with much better time resolution than before 
(future monitoring observations are desirable). Detailed studies of the mass transfer 
process (stellar wind, magnetosphere) can be performed as well (e.g., 4U 1907+09, 
Vela X-l, Cen X-3). Future work is foreseen to include a comparative study of the 
Suzaku sample applying the new cyclomc line model. 

Acknowledgments: We thank the Suzaku team for their help in executing these observations, espe- 
cially for the excellent timing of the 1A 1118— 61 ToO (Fig. 1, left) and the eclipse-to-eclipse Cen X-3 
observation (not discussed here due to space limitations). 
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